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Energy spectra of a particle with mass m and charge e in the cubic Aharonov-Bohm billiard con- 
taining around 10 4 consecutive levels starting from the ground state have been analysed. The cubic 
Aharonov-Bohm billiard is a plane billiard defined by the cubic conformal mapping of the unit disc 
pervaded by a point magnetic flux through the origin perpendicular to the plane of the billiard. 
The magnetic flux does not influence the classical dynamics, but breaks the antiunitary symmetry 
in the system, which affects the statistics of energy levels. By varying the shape parameter A the 
classical dynamics goes from integrable (A = 0) to fully chaotic (A = 0.2; Africa billiard). The level 
spacing distribution P(S) and the number variance S 2 (L) have been studied for 13 different shape 
parameters on the interval (0 < A < 0.2). GUE statistics has proven correct for completely chaotic 
case, while in the mixed regime the fractional power law level repulsion has been observed. The 
exponent of the level repulsion has been analysed and is found to change smoothly from to 2 as the 
dynamics goes from integrable to ergodic. Further on, the semiclassical Berry-Robnik theory has 
been examined. We argue that the semiclassical regime has not been reached and give an estimate 
for the number of energy levels required for the Berry-Robnik statistics to apply. 

PACS numbers: 05.45.+B, 05.40.+J, 03.65GE, 03.65.-W 

The objective of this paper is to study the energy level 
statistics of a system without antiunitary symmetry. In- 
troducing magnetic field to a plane billiard violates the 
time reversal symmetry. However, if there are space sym- 
metries present, there can be a combined space-time sym- 
metry in the system, which is antiunitary [Q. Therefore, 
a billiard with no space symmetries is required. The nat- 
ural choice is the cubic billiard defined by the conformal 
mapping of the unit disc (\z\ < 1) 

w = z + Xz 2 + Xe^z 3 (1) 

onto the physical plane w, with the shape parameter A. 
A point magnetic flux through the origin perpendicular 
to the plane of the billiard is introduced. Such billiard 
systems are called Aharonov-Bohm billiards. The clas- 
sical dynamics is unaffected by the point fluxQ and the 
system has the scaling property, which says that the clas- 
sical dynamics is the same at all energies. The particle is 
therefore considered to have a unit speed. This family of 
billiards was introduced and first studied by Berry and 
Robnik §]. 



*Only the trajectories which pass through the origin ex- 
actly are affected, but they have measure zero among the 
trajectories. 



Energy level statistics in the limiting, fully chaotic case 
of Africa billiard (A = 0.2), has been found to be con- 
sistent with the Gaussian Unitary Ensembles (GUE) of 
random matrices . In the present work we have inves- 
tigated the energy level statistics in the regime of mixed 
classical dynamics and have found the fractional power 
law level repulsion]^, while we argue that the semiclassical 
theory of Berry and Robnik Q applies for much larger 
sequential quantum numbers (see also Jl6{), in our case 
say 10 8 . We also reconfirm [^|] the quadratic level repul- 
sion and the validity of GUE (level spacing distribution 
and number variance) statistics in the ergodic case with 
high statistical accuracy. 

The Schrodinger equation for the system ^ is 



2 



r z o(f> r z 



dw 



dz 



2 



ip{r,<j>) = , 



^(1,0 = 0. (2) 
Here k 2 = 2l jp- where E is the energy and the parameter 



^The phenomenon of fractional power law level repulsion has 
been typically found in mixed systems with antiunitary sym- 
metry and not too high energies (see and the references 
therein) . 
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a = where $ is the magnitude of the magnetic flux 
and 2irh is the Planck's constant. By (r, </>) we denote the 
polar coordinates in the z plane. The relevant region of 
the values of a is [0, |]. It can easily be shown |l]-|3| that 
the antiunitary symmetry is still present for integer and 
half integer values of a, so the most appropriate value is 
a = j. We choose exclusively this value in our present 
calculations. In the integrable case (A = 0) we expect to 
find Poissonian energy level distribution 



P(S) 



(3) 



while in the completely chaotic case P(S) is predicted 
by the GUE of the random matrix theory. The Wigner 
approximation ^\ gives for GUE 



P(S) 



(4) 



For small S we have P(S) oc S 2 , i.e. quadratic level 
repulsion. 

Let us first describe a semiclassical argument (see |J) 
to estimate the number of converged energy levels in our 
calculation. By "converged" we mean the accuracy of 
at least one percent of the mean level spacing. We ex- 
pressed the matrix elements of the Hamiltonian H in a 
basis defined by eigenvectors of the integrable Hamilto- 
nian H° (circular billiard, A = 0, see [p|]lC|]). We trun- 
cated the basis at N°th vector and diagonalised the fi- 
nite (N° x N°) matrix. The eigenvectors |n) of H° are 
defined as H°\n) — e n \n) and the eigenvalues of H as 
H\a) — E a \a). The question is, how many accurate en- 
ergy levels E a associated with H are there in such calcu- 
lation? In the semiclassical limit the Wigner transform of 
an eigenstate with energy E is localized on the classical 
energy surface (see and the references therein). 




FIG. 1. Number of converged levels as a function of the 
shape parameter A. 

This yields that an energy level E a of H will con- 
verge if all eigenvectors \n) of H° with the energy sur- 
faces H°(q,p) = E® intersecting the energy surface of H, 
H(q,p) = E a , are present in the truncated basis. [] Geo- 
metrical consideration of foliation of energy surfaces (see 
and the references therein) gives 



N 
TV" 



r(A) 



A(X) 



C|*|<1 



I dw | 
I dz I 



(5) 



where A(X) is the area enclosed by the billiard bound- 
aries and is equal to A(X) = ir(l + 5A 2 ) for the cubic 
billiard. The fraction of the converged levels as a func- 
tion of A is plotted in Fig. [l| Approximately one third 
of the levels are correct at A = 0.2, whereas at A = 0.05 
there are about two thirds of good levels. As mentioned 
above, by a good level we mean the accuracy of at least 
one percent of the mean level spacing. Of course, this a 
priori criterion is helpful, but the final judgement of the 
precision was based on the actual numerical convergence 
of the levels. 

If the basis functions are arranged in a smart way, the 
jyo x jy-o ma trix which has to be diagonalised has band 
structure. In fact we do not diagonalise the Hamiltonian, 
but its inverse H~ l Q which can have the desired band 
structure. The rearrangement of the basis is described 
in ||. The CPU time demand of the diagonalisation 
drops from 0(N 3 ) to 0(N 25 ) due to the rearrangement. 
The size of the matrix was N° = 18000 for A = 0.2 and 
N° = 15500 for thirteen other shape parameters from the 
Table |[ By diagonalising the band matrices the fourteen 
spectra have been obtained containing r(X)N° accurate 
levels. 

Using the Weyl formula the levels have been unfolded 
to unit mean level spacing and nearest neighbour level 
spacing distribution P(S) and number variance E 2 have 
been calculated. Instead of P(S) the cumulative level 
spacing distribution W(S) = J Q P(s)ds is shown in Fig. 
||. There are the so-called U -functions shown in the 
insets. By fitting data locally (S < 0.5) to the Brody dis- 
tribution, pP(S) = aS 13 exp[~bS l3+1 ) @, the behaviour 
of P(S) for small S in the transition region < A < 0.2 
has been analysed. However, we know that the Brody 
distribution can not globally (for all S) capture the phys- 
ical P(S) (at least for > 1). Nevertheless, we found a 
good agreement with the Brody distribution for S < 0.5. 



*We tested the quality of the Wigner approximation by com- 
paring it to the exact infinitely dimensional solution (see 
and found that Wigner approximation is almost perfectly cor- 
rect, the tiny deviations from the exact solution are two orders 
of magnitude smaller than the deviations of our numerical 
results. 



nt can be shown that the terms in the Hamiltonian (|2|) con- 
taining a are of the second order in h, so the magnetic flux 
does not influence the density of states in the semiclassical 
limit (of course, it mixes the levels and affects the statistics) . 
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FIG. 2. The cumulative level spacing distribution W(S). 
Numerical data: full line, Poisson: dotted line, GUE: dashed 
line and Berry-Robnik: dash-dotted line. For the definition 
of the U functions plotted in the insets, see text. 




FIG. 3. The function T(ln S) for some of the calculated 
spectra. The value of the shape parameter is written in the 
upper left corner of each figure. 

The results are shown in Fig. || and in Table | and 
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Fig. ||. In the Fig. || there is the T function shown, 
which transforms Brody distribution to the straight line: 
TQnS) = ln[-ln{l - W(exp(lnS))}]. Fitting the nu- 
merical data to the line, the slope is the exponent of the 
level repulsion plus one. The statistical error of the mea- 
sured values is shown with the bars. 

A good agreement with the Poisson statistics is seen in 
the integrable case (A = 0), while small deviations from 
GUE statistics in the fully chaotic case are observed. The 
deviations at A = 0.2 are probably due to localization 
of the wave functions, e.g. scars, and further work is 
necessary to explain that. For small S the quadratic level 
repulsion is verified at A = 0.2. The fractional power law 
level repulsion with the exponent /3 varying continuously 
from to 2 is observed. From the Fig. y we can also 
confirm that the Brody distribution is valid only locally, 
for S < 1. 

As for the number variance E 2 (L), we have observed 
continuous transition from Poisson, Y? Poisson (L) = L, to 
GUE, Y? GUE {L) w 0.10 hxL + 0.34, statistics. The satu- 
ration sets in at L* ~ 20, which is roughly in agreement 
with Berry's semiclassical theory . In the fully chaotic 
case we have found very nice agreement with the GUE. 

We also explored the relevance of the Berry-Robnik 
theory, which is based on the statistically independent 
superposition of Poissonian sequence of regular levels and 
GUE sequence of chaotic levels, each having statistical 
weight equal to the relative measure of the corresponding 
classical invariant component. 

Following the procedure in JjJ we derived the Berry- 
Robnik formula for Poisson-GUE transition: 



P ati (S) 



-PiS 



+ (2pip 2 - p\p2S) erfc(-^p 2 S) 



+ 
(6) 



E(L) 



where p\ = 1 — p 2 is the relative measure of the regular 
part of the classical phase space. The classical work is 
to determine pf and pf- Different methods and difficul- 
ties connected to this work are presented in p!Bj . After 
the careful study the following method has been chosen. 
The SOS is first divided into a mesh of M x M cells. 
A trajectory is then run on the largest chaotic compo- 
nent. Each cell i has a counter which counts the number 
of times a trajectory passes the cell, rjj. After having 
made enough steps the distribution of the numbers r\ is 
plotted. The distribution has a peak around i) = 0, a 
minimum and another maximum for larger values of r\. 
If the distribution is normalised, the area to the right of 
the minimum is equal to p 2 . In table 1 there are values 
of pi for different shape parameters A. 




FIG. 4. The number variance statistics S 2 (L) for different 
shape parameters. The limit cases are shown dotted (Poisson) 
and dashed (GUE) 



A 


cl 
P2 


BR 
P2 


P 


0.000 


0.00 


0.18 


0.00 


0.010 


0.05 


0.26 


0.02 


0.017 


0.15 


0.28 


0.03 


0.022 


0.26 


0.32 


0.05 


0.025 


0.34 


0.35 


0.09 


0.028 


0.58 


0.36 


0.03 


0.033 


0.75 


0.41 


0.15 


0.040 


0.86 


0.49 


0.16 


0.049 


0.92 


0.68 


0.48 


0.055 


0.99 


0.78 


0.69 


0.062 


0.99 


0.77 


0.87 


0.070 


0.99 


0.94 


1.42 


0.100 


1.00 


0.97 


1.58 


0.200 


1.00 


0.99 


1.93 



TABLE I. The classical value pf for different shape param- 
eters A. The third and the fourth column are the best-fitting 
values for the Berry-Robnik and Brody model, respectively 
(see further text). 



"'Surface of Section, chosen in the standard way, see Q and 
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FIG. 5. Calculated values of pf, p2 R an d ^PBrody as func- 
tions of A. 

The quantum parameter p BR is determined through 
fitting to the equation ra) and should be compared to 
the classical measure p% , The results are presented in 
Fig. U and in Table |. In insets of Fig. [| the fine 
scale deviations from the Berry-Robnik formula are plot- 
ted in terms of the U(W(S)) - U(W BR {S)) vs. W(S) 
(full line). We also plot U(W Brody (S)) - U(W BR (S)) for 
comparison (dashed line). The transformation U(W) = 
— arccos \fl — W is used pj in order to have uniform sta- 



tistical error over the plot, SU = l/(n\/N). The dotted 
horizontal lines indicate ±SU. The Berry-Robnik fit is 
not statistically significant and the values of the fitting 
parameters p^ R are different from the classical values pf. 
Why this is so and how the crossover from Brody- like be- 
haviour to Berry-Robnik takes place is discussed in [jl6| . 
As one can see at small S (and therefore at small W(S)) 
the agreement with Brody is better than Berry-Robnik, 
and occasionally, like e.g. at A = 0.055 and 0.040, even 
globally Brody seems to capture the behaviour of our 
data, although, on the other hand, there are expected 
deviations at large S (and large W(S)) e.g. in the fully 
chaotic case of Africa billiard (A = 0.2), where (3 — 1.93, 
and Brody asymptotic behaviour at large S is certainly 
wrong, because its exponent varies oc S 3 , which is dif- 
ferent from GUE in equation (^), oc S 2 . What we see 
in this [/-function representation is consistent with the 
T-function representation in Fig. || 

In this paper we have presented the numerical computa- 
tion and statistical analysis of energy spectra of a Hamil- 
tonian system with mixed classical dynamics and with- 
out antiunitary symmetry, namely the Aharonov-Bohm 
cubic conformal plane billiard. We reconfirmed the GUE 
statistics in the ergodic (fully chaotic) case with approx- 
imately 6000 good consecutive energy levels whose accu- 
racy is better than one percent of the mean level spacing. 
Further, in the KAM regime we have found the fractional 
power law level repulsion with the exponent varying 
smoothly from zero to two as the corresponding classical 



dynamics goes from integrable to ergodic. On the other 
hand we know that in the ultimate semiclassical limit the 
Berry-Robnik theory |?],[l6| should apply, which does not 
exhibit level repulsion (P BR (S = 0) = p\ ^ 0). The 
resolution of this disagreement lies in the fact that the 
key assumptions of the Berry-Robnik theory arc fulfilled 
only for extremely large sequential quantum numbers or 
small effective Planck's constant. Namely, one should re- 
quire extendedness of the Wigner functions on the whole 
corresponding classically allowed invariant components 
of phase space. This is true only if the Heisenberg time, 
tji = 2irh/ (AE) (the time until the discreteness of the 
spectrum is not yet resolved and the quantum evolution 
follows the classical one), is larger than the classical er- 
godic time, t e (the typical time after which classical dis- 
tributions reach equilibrium steady state on the invariant 
set). In our case we have typically t e > 10 4 (velocity is 
one), so the condition tn > t e yields that the sequential 
number N should be larger than about 10 8 . We believe 
that our present work contributes to the new develope- 
ments in the wide field of research of quantum chaos, 
recently reviewed by Weidenmiiller et al. p7[. 
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